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Abstract 


Although heavy ion transport theory is developed to a relatively 
advanced stage, the present limitation in biomedical and electronic 
applications is the uncertainty in nuclear fragmentation parameters. 

Present status on 20 Ne beams is discussed and useful formulae are presented 
for future use in analysis of beam transport experiments. 
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Introduction 

Considerable progress has been made, in recent years, toward the 
understanding of the interaction of high charge and energy (HZE) particles 
with extended matter. A number of theoretical methods have emerged by 
which estimates of transport quantities can be made from microscopic 
interaction data . 1 " 3 Although the microscopic interaction data have been 
poorly known, improvements in quantity and quality are becoming available 

o c 

so that meaningful calculations could be made in the near future. ~ 

In the present paper, we develop a simplified theory of HZE transport 
which is suitable for analysis of ions of modest energies. The main value 
of the formalism is that a simple framework is provided by which all 
important transport quantities can be calculated from the basic 
interaction data using a small desktop computer. We will show that major 
uncertainties in transport quantities arise due to uncertainties in 
present-day fragmentation parameters. The present theory may be useful in 
preliminary data analysis in addition to its importance as a step in 
developing more detailed and accurate transport models. In the present 
analyses, 20 Ne beams are investigated due to the availability of 
experimental data. 

Transport Theory 

With a straightahead approximation, the basic equations are given by 3 

li-li ~V E)t< ’j]*j (x ’ E > ’ ^kVk'*^ O) 

where the nuclear cross sections are assumed energy independent and 
approximately given by 7 

crj = 58.1(A 1/3 +2.51-l.ll/A 1/3 -l.ll/2.51) 2 


( 2 ) 
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The stopping power S j ( E ) is estimated using the approximate range energy 
relation given by 

Rj(E) = «E n yZj 2 (3) 


where 6 = 9.2 x 10 -4 , n = 1.75, Aj is the ion mass, Zj is the ion charge 
and R j ( E ) is the ion range in cm of G5 nuclear emulsion for an assumed 
emulsion density of 3.815g/cm 3 . For a monoenergetic beam of I-type ions of 
energy E 0 incident on the emulsion stack the HZE fluence is given by 8 


S.[R.(E)+x] . 

Mx,E) = e -a j x <j>j {0,Rj _1 [Rj (E)+x]} 


V E) 


R 1 Cx+R 1 (E)3 exp[-o.R. ( E " ) ] 

+ -=— fJ J Jm. |/ a l/ ^[x+R.(E)-R.(E ,, ),E"]dE" 

S .(E) E exp[-ajRj(E)] Z J k k * kL ^ J V 

J 


(4) 


where 

+j {°» R j“ 1 C R J ( E) + x 3} - « { R j " ! CR j { e ) +x3- e q } i (5) 


which are solved by an iterative procedure. The first term, found by neglecting 
the integral term, is 



(x , E ) = 


e-V 

V E) 


6 iI 6[x+R j (E)-R j (E 0 )] 


( 6 ) 



4 


The second term is derived by substituting 4>j (°) (x,E) and evaluating the 
integral in equation (4) as 

m m,,o. R.^CMEJ+x] 

. (^(x.E) .-iL-L / J J exp{-0.x -0 I [ME)-R I (E”)]-a i [R.(E-)-ME)]( 

j S.(E) E 1 1 J J J 

J 

d[x + R j 4E)-R j (E") + R I (E M )-R I (E o )]^ FT (7) 


This is evaluated to yield 

<K- (1) (x,E) = expf-ajX-ajCRJEJ-R (E")] 

J S J< E) Ivi-vjI 

- °j C R j (E")”Rj (E)]} 

where vj is the ratio Zj 2 /Aj and Rj(E") given by 
R,(E") = [R r (E )-R . (E)-xj 

J ' ' V * V • I ' o' J v ' J 

J 

is subject to the restriction 

Rj(E) < Rj (E") < x +R j ( E ) . 

Note that equations (9) and (10) result in 

^■[Rl(E 0 )-x] < Rj (E) <^R,(E 0 ),x 
j J 


(R) 


(9) 


( 10 ) 


t 


(11) 
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We now consider evaluation of the dose, the integral fluence, and the integral 
LET spectra. 

The dose of the primary beam is 


D I ( 0 ) (x) = A f 4,/ 0 ) (x,E) S (E)dE 
o 

- A,Sj (R ] - 1 [R I (E 0 )-x])exp(-o,x) (12) 

and similarly 8 

D j (1)(x) = A j f V E),|, j {1)(x ’ E)dE 

v.m. T a T “ 

= A. / expf-ajX-ajrR^Fj-R^E 11 )] 

rr v il ° 

- o j [R j (E")-R j (E)]}dE (13) 


where the relations (9) and (10) are yet to be employed. The limits of 
integration are to be restricted by equation (11). The integrand of 
equation (13) has the values 


exp (-o jx) 
exp(-a.x) 


(E" 


(E" 


E) 


} 

_1 [x+Rj ( E ) ] ) 


(14) 


Approximating the integrand of equation (13) with an exponent which varies 
linearly between the limits indicated by equation (14) yields 



m ji g iVj 

v r v. 



( CT r a j> x 


[exp(-a j x)-exp(-a I x)] 


(15) 
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as given elsewhere. 9 In (15), E». and E u . are the limits of E given by 
equation (11). 

The integral fluence is similarly derived as 


with 


♦ , (0) (x) = ; ^""(x.EJdE - e-I 


(o), 


-cr T x 


*j ( 1 ) (x) 


R , , 

v.m, T a T u i (cr--a T )x 

= 77 TT l ex P{~ CT i x ~ r .ft ( R j- R . )}d R 


V v j! R £ 


U j 


(16) 


R -R 

= V j m j lC?I 1 U J £ j ^ Cexp(-a.x)-exp(-OjX)] 


I I (a T -a".)x 

rr v jl 1 J 

and where we have made use of equation (14) and the values 


v r 

R . =~r R r(E n )-xl 
i. v. L I v o' 

J J 


R = _L R T (E )-x 
U . v. I' o' 

3 3 


} 


(17) 


(18) 


given as the limits of equation (11). 

The LET spectrum is evaluated in a manner similar to equation (17) as 


*j(x,>S) 


= / <f>.(x,E)dE 


(19) 


m 


where E$ is the energy at which Sj(E) is numerically equal to S and is 
the energy at which Sj ( E ) obtains its maximum value. In terms of equation (3) 


E$ - [AjVj/n6S] 


1 

.n=T 


(20) 
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Defining R$, = Rj(E s ) and (L.= Rj (E m ) we may write equation (19) 
as 


/ i \ v .m . T a T 

*, (1) (x,>S) --ULI 

J V 


R -R „ 
u . i. 


ii-i • -i v (^T-a-)x 

— - r ■ - - "-4 ? {exp[-a.x - j* l—, (R -R- )] 

r vj| {a r a j )x J ( u j m j j j 


- expC-CTjX) } 


( 21 ) 


The maximum stopping power in the present calculation is taken as 

S.(Em) = 10 3 Z 1,1 (22) 

J J 

as determined from the results of Zeigler. 10 

Results 

The transition of a 20 Ne ion beam of 250 MeV/amu in G5 nuclear emulsion has 
been studied experimentally by Jain and Agyarwal. 7 At each nuclear 
fragmentation site, the number of fragments of each species were identified by 
observing the 6-ray density along the track. A set of fragmentation cross 
sections was compiled. The calculated dose given by equation (15) using Jain's 
and Aggarwal 's fragmentation cross sections are shown in figure 1 in comparison 
to their experimental values of dose from various secondary components. As can 
be seen in figure 1, the dose obtained from counting particle tracks in emulsion 
is reasonably consistent with their fragmentation parameters. Additional 
fragmentation cross sections are shown in table I. The values labeled TS are 
the Tsao-Si 1 berberg parameters calculated using their computer programs. 11 The 
values labeled VR are the velocity renormalized Tsao-Si 1 berberg values which 



were found to give the best agreement with 20 Ne transition experiments 
performed by Schimmerling and coworkers 3 at the Lawrence Berkeley 
Laboratory. The values labeled LaRC are recent results from theoretical 
calculations at the Langley Research Center. 6 

In figure 2, the Jain and Aggarwal dose curves are compared to 
evaluation of equation (15) using the fragmentation parameters of Tsao and 
Silberberg. It is seen that the N,0,F group greatly overestimates the 
corresponding experimental values. The Li,Be,B,C group is only slightly 
underestimated by the Tsao-Si Iberberg parameters. The He dose is 
underestimated by both sets of parameters. It is clear from table I that 
the VR parameters would make an even larger overestimate of the 
experimental values for the important NOF group. 

In an earlier study of 670 MeV/amu 20 Ne ions in water, we found that 
the Tsao-Si Iberberg parameters significantly underestimate the dose at a 
given depth. 3 The velocity renormalized parameters (VR), however, provide 
good agreement with the LBL experiments as seen in figure 4 of the earlier 
study. It is clear from comparison of results for the Tsao-Si 1 berberg 
fragmentation parameters with results for the Jai n-Aggarwal parameters that 
if used to perform the Ne on water calculations the latter values would 
lead to even further reductions in dose and thus lead to even greater 
disagreement with the LBL experiments. 3 This is especially clear since the 
velocity renormalized values in table I are in good agreement with the LBL 
experiments . 

The present situation is most unclear with respect to the 
fragmentation and transport of 250 MeV/amu 20 Ne ion beams in extended 
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materials. The results of Jain and Aggarwal are internally consistent but 
cannot be reconciled with the experiments of Schimmerling et al. 3 insofar 
as the transport equation (1) is correct. It does not appear that the 
simple range formula is responsible for these differences since more 
accurate values were used in our earlier results. 3 An additional set of 
experiments on 20 Ne at 250 MeV/amu using sophisticated experiments such as 
those of Schimmerling et al . is likely to settle the issue. 

The calculated LET spectra is shown in figure 3 for the three 
fragmentation data sets of Jain and Aggarwal, Tsao and Silberberg, and 
velocity renormalized. Only the high LET components are shown in figure 3 
which demonstrate an order-of -magnitude uncertainty in calculated LET 
spectra due to differences in fragmentation cross sections. Such large 
uncertainty in LET spectra will not be acceptable in the evaluation of 
health hazard/benefit as well as estimation of effects on active 
electronics devices. Clearly, the principal weakness in developing 
transport models is the lack of consistent fragmentation parameters. 

To help define adequate fragmentation cross sections, the Langley 
Research Center is supporting a theoretical program to evaluate the 
fragmentation parameters from fundamental nuclear theory. 6 Preliminary 
results of these calculations are shown in table I for comparison. It is 
hoped that such physically based models, coupled with the relatively 
limited experimental data base, will one day allow reliable predictive 
capability in heavy ion radiation research. 

It is hoped that the present formulas will also be helpful in 
understanding the transport of heavy ions through extended matter. Higher 



order terms can be added to the series approximation when necessary. The 
greatest uncertainty, however, in most applications will lie in the 
inherent uncertainties of the fragmentation parameters rather than in the 
basic transport formalism presented. 
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Figure 1.- Experimental dose of Jain and Aggarwal compared 
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